Exploiting structured environments for efficient energy transfer: 
The phonon antenna mechanism 
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A non-trivial interplay between quantum coherence and dissipative environment-driven dynamics is becoming 
increasingly recognised as key for efficient energy transport in photosynthetic pigment-protein complexes, and 
converting these biologically-inspired insights into a set of design principles that can be implemented in artificial 
light-harvesting systems has become an active research field. Here we identify a specific design principle - 
the phonon antenna - that demonstrates how inter-pigment coherence is able to modify and optimize the way 
that excitations spectrally sample their local environmental fluctuations. We place this principle into a broader 
context and furthermore we provide evidence that the Fenna-Matthews-Olson complex of green sulphur bacteria 
has an excitonic structure that is close to such an optimal operating point, and suggest that this general design 
principle might well be exploited in other biomolecular systems. 
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Introduction. — Experimental techniques such as optical 2D 
Fourier Transform spectroscopy have recently begun to probe 
the ultrafast photophysics of energy transport in a range of 
pigment-protein complexes (PPCs) taken from photosynthetic 
green sulphur bacteria, marine algae and higher plants EH?). 
All of the key photosynthetic light reactions (photon capture, 
energy transport and charge generation) are carried out in PPC 
structures |5 , 6], and the often > 90% quantum efficiency with 
which they carry out these functions has created considerable 
interest in understanding and exploiting their design features 
for artifical solar energy applications 0171. Surprisingly, re- 
cent experiments have found direct evidence for long-lasting 
quantum coherence amongst the excitons which transport en- 
ergy in these complexes. In the case of the Fenna-Matthews- 
Olson (FMO) complex, these coherences can persist on pi- 
cosecond timescales in cryogenic conditions and are still ob- 
servable at room temperatures 0. Several phenomenolog- 
ical theories have subsequently shown that there is an opti- 
mal mixture of coherent inter-pigment energy transport and 
stochastic environmental noise that may lead to faster and 
higher-yield energy delivery in PPC architectures, suggesting 
that quantum effects may underpin their efficient function (9]- 
[T5i Microscopic investigations have also recently shown the 
key role of both discrete and continuous environmental fluc- 
tuation spectra in stabilising the long-lasting coherences ob- 
served in spectroscopy as well as facilitating efficient trans- 
port fT6ti2Ql . 

In this letter we build upon these microscopic studies to 
describe an explicit design principle through which coher- 
ent excitonic coupling is exploited to optimize energy trans- 
port in the typical environmental conditions found in PPCs. 
This is based on the idea that coherent couplings can allow a 
pigment network to be spectrally 'tuned' into configurations 



where the emerging relaxation (transport) pathways may ex- 
tract maximum noise strength from the environment and thus 
proceed faster. In analogy to the well-established idea that 
excitonic interactions in light-harvesting antenna complexes 
allow pigments to absorb a broader range of the ambient light 
spectrum |6l, we refer to this principle of advantageous, 
coherently-modified sampling of the protein noise spectrum 
as the phonon antenna mechanism fT9l . 

The phonon antenna. — To illustrate the concept, we begin 
by setting out the (canonical) model of exciton transport in 
a PPC network whose function is to deliver excitonic energy 
to a reaction centre (where it is used to generate electrons for 
subsequent - dark - stages of the photosynthesis) (6). The net- 
work consists of a set of chromophores (sites) which are held 
in space by a protein matrix. Each chromophore possesses 
an optically excited state | i) , whose energy £ ; depends on the 
local environment. The electronic Hamiltonian Hg is then 
given by H s = |/) (/| + Li^jJij 10 01. the sum of single 
exciton contributions on each site i with local energy energy 
£i and coherent coupling terms describing the inter-site ex- 
change of excitations with dipolar coupling strength Jij. The 
states obeying Hs\e n ) = E n \e n ) are in general delocalized ex- 
citonic states, which can be also be expressed in the site basis 
as \e n ) = L/QIO- The ultimate transfer of the exciton to a re- 
action center is modelled by including an isolated site (labeled 
and referred to as the 'sink') which is populated irreversibly 
from a particular local site l9lfTTll. 

Non-electronic degrees of freedom, due to the surround- 
ing protein environment and/or internal vibrational modes of 
the chromophores, also couple to the excitons, inducing lo- 
cal fluctuations of the site energies These 'environmen- 
tal' fluctuations (phonons) are modelled as originating from a 
continuous bath of harmonic oscillators coupled to the exci- 
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FIG. 1. (a) A model system consisting of a three-site network show- 
ing coherent couplings (double-headed arrows) and relative local site 
energies of state |/) subject to dephasing noise as described in the 
text. An excitation on site 1 will be transferred to the sink (site 0) via 
sites 2 and 3. The population transfer from 3 to the sink is incoherent 
(dotted arrow), (b) Diagonalising the network Hamiltonian results 
in new eigenstates \ei). For the parameters in (b), the energy differ- 
ences between eigenstates correspond to maxima of a hypothetical 
environmental spectral function shown in (c). Environment-induced 
relaxation between these states (single-headed arrows) is thus opti- 
mised, as transition rates are proportional to the spectral density at 
the exciton energy differences. 

ton and constitute the primary source of energy relaxation and 
dephasing for the excitonic system (6l [22) - see Appendix A. 
Coherent inter-site couplings in Hg can be stronger or of the 
same order than the characteristic strength of the environmen- 
tal interactions |Q]|3llU|23). m suc ^ cases delocalised exciton 
eigenstates, for which Hg is diagonal, therefore provide the 
preferred basis to analyze the system's transport dynamics. 
Within this framework, we can employ secular Redfield the- 
ory to obtain the following Markovian equation of motion for 
the population p nn of the exciton states {\e n )} I24ll25l : 

Pmm W mn p nn ^ ( W nm ~\~ ^e m ) Pmm •> ffl^O 

POO = ^en^Opnn, (1) 

where the transition rates between eigenstates are given by: 

= f 27Z ( /((D mn )% mn n((D mn ), CO mn > 

mn I 2n f ((D nm )% mn [n((D nm ) + 1], GW<0 

with 

Xmn ~ I C n C m | , n((D) — i" 5 (Djnn = E m E n . 

i e k B T 

Here the irreversible process that locally populates the sink 
from site n is assumed to have an (environment-independent) 
rate = 1/(1 P s ) EHIIL which is typical for transfer to 



reaction centers. In the exciton picture this leads to popula- 
tion transfer to the sink at a rate T en ^$ = | (e n \n) | 2 r^o- Our 
figure of merit for assessing the transport efficiency of a given 
network is the population in the sink p S u±(t) = Poo(0 at some 
time t following the local excitation of a selected site. The 
population dynamics is entirely driven by the environment 
and, as a result of energy conservation, the relaxation rate of 
an exciton transition \e n ) —> \e m ) (which requires the emission 
of a phonon) is thus dependent on both the phonon density 
of states and the coupling strength to phonons at frequencies 
equal to the energy difference E nm . Both these properties are 
described by the environmental spectral function ^ ((D) and 
the transition rate W nm is proportional to J? (E nm ). As typical 
PPC spectral functions often have (at least) one maximum at 
a finite frequency ifTTI and may exhibit sharp features due to 
long-lived vibrational modes in the environment, it becomes 
possible to tune the Hamiltonian parameters to create exci- 
tonic states with energy differences E nm at local maxima of 
the spectral function and thus optimise W nm . This is the basis 
of the phonon antenna effect. 

Figure [T] illustrates the phonon antenna mechanism in a se- 
quential transport scenario. The simplest model consists of 
a network made up of three sites, |/), with local energies 
z = 1, 2,3, as shown in Figure[T^. Sites 1 and 2 are assumed to 
be nearly degenerate and strongly coupled, with a smaller cou- 
pling between sites 2 and 3, and almost no coupling between 
sites 1 and 3; that is, /13 = 0,/23 < J\2 < £1 ~ £2 = £, £3 = 0. 
Initially, only site 1 1 ) is excited. For these parameters, the 
coherent coupling J12 splits the degeneracy of states |1) and 
1 2) and the exciton eigenstate of Hg are approximately \e\) w 
2- 1 /2(|i) + | 2 », \e 2 ) « 2-^(1 1) - |2» ,h> « |3>, with en- 
ergies E\ — £ -\- J 12^2 = £ — hi and £3 = 0. The initial con- 
dition can be re-expressed as |1) = 2~ l ^ 2 (\e\) -\-\e2}), and an 
energy diagram for these new states and initial condition is 
shown in Figure [TJ). The efficiency of transport to site 3 (and 
thence to the sink) now depends on the dissipative rates W\2 
and W23 of the environment-induced transitions \e\) — » ^2) 
and ^2) — > 1 3), respectively. 

Let us consider now an spectral function possessing a max- 
imum at an energy co#, as depicted in Fig. [I]:). When 
J\2 = £/2 = (Oh, we find that £"12 — E23 ~ (Or and a ladder 
of equally spaced states appears with a splitting (Oh- Thus 
a clever choice of J12 simultaneously optimizes both transi- 
tion rates W\2 and W23 and generates efficient energy trans- 
port to the sink. This simple system illustrates the essen- 
tial idea behind implementing an efficient phonon antenna: 
Electronic interactions are tuned (e.g. by controlled or evolu- 
tionary adaption of distances and relative orientation ofchro- 
mophores) to achieve resonance with maxima of the environ- 
mental vibrational spectra. 

This already suggests how in more general networks, co- 
herent couplings and site energies could be tuned to enhance 
(or inhibit) a particular set of energetic transitions (which in a 
multi-chromophore network could also correspond to a partic- 
ular spatial pathway). The more structure in the spectral func- 
tion, the greater the possibilities for phonon antenna effects 
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to be harnessed, as shown in Fig[T]: where a shoulder feature 
in ^f(co) also enhances the direct W\ 3 transition. In anticipa- 
tion of our analysis of the FMO complex, we now perform a 
quantitative analysis of this three-site model system. 

Motivated by actual values in typical PPCs, the consid- 
ered Hamiltonian parameters are, in units of cm -1 : £\ = £2 = 
300, J l2 = J21 = 100, J 23 = J32 = 30, £ 3 = J 13 = J 31 = 0. The 
spectral density function is taken to be of a quasi-Lorentzian 
shape centered in (Oh (26) : 



e 2 (cm *) 



CO) 



{(0 2 -co 2 ) 2 + {2nTco) 2 



(3) 



where j3 is fixed in order to have a reorganization energy X c± 

K(0 "^ = 35 cm -1 , and T = 60 cm -1 . We assume that initially 
site 1 is selectively excited. 

Within this parameter regime, the inset in Figure [2] shows 
the population in the sink after 2 ps as (Oh is varied. As 
expected, a clear peak in the transport efficiency emerges 
when the spectral function is peaked slightly above the value 
(Oh = 200cm -1 which would correspond to perfect match- 
ing of the inter-exciton energy differences to the maximum 
of the spectral density. Figure [2] shows the sink populations 
after 2 ps for fixed (Oh and £\, while varying the coupling 
J12 and the site energy £2. We observe a symmetric structure 
in the transport efficiency with respect to J 12 —> —J11, but 
find that the degenerate configuration £\ = £2 is not a global 
maximum of the efficiency. The highest efficiency occurs for 
£2 ~ 240cm -1 and J\2 ~ 95cm -1 . This larger value of the 
efficiency is the result of a slightly enhanced coupling (over- 
lap) between sites 2 and 3 that results from shifting site 2 to 
a lower (local) energy. However, the differences in exciton 
energies E\—E2 = 200cm -1 and E2—E3 = 180cm -1 remain 
close to the maximum of the spectral density. Indeed we no- 
tice that the reduction in the optimal values of J\2 as £2 de- 
creases maintains energy differences close to the spectral den- 
sity maximum. The shift in the global maximum illustrates 
the role of the % factor in determining the energy and cou- 
pling structures which optimally exploit the functional form 
of the environment's spectral function. 

The phonon antenna in the FMO complex: Modified Red- 
field Simulation. — The FMO complex is a PPC present in 
green sulfur bacteria which acts as a 'molecular' wire by me- 
diating exciton transfer from the light-harvesting chlorosomes 
(antenna complex) to the bacterial reaction center (5l[6). It has 
a trimeric structure with each of the three monomers contain- 
ing seven chromophores bound to a common protein struc- 
ture. Under low light conditions, it can be assumed that there 
is only one excitation at a time, and that all monomers oper- 
ate independently. The excitation is thought to originate at an 
eighth site (common to several monomers) which injects an 
exciton locally at sites 1 or 6 of a given monomer [27]. The 
exciton then relaxes through the manifold of exciton states, 
and is directed towards site 3 (which is the primary compo- 
nent of the lowest energy exciton state). Due to its proximity 
to the reaction centre (RC), the exciton on site 3 is then rapidly 
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FIG. 2. Transferred sink population after 2 ps for the Markovian 3 
site model system described in the text and depicted in Figure 1. The 
environment is described in terms of a Lorentzian spectral density 
peaked at a frequency &># = 200cm -1 for a temperature T = 4K. 
Actual values of the Hamiltonian values, lie as expected (blue dot), 
at the top of a maximum in transfer yield. In the inset, population 
in the sink after 2 ps for a spectral density peaked at given a 
temperature T = 4K. As an example were the matching takes place, 
the phonon antenna is formed and adjusted to a peak &># ~ 2J = 
200 cm -1 . The other peak around (Oh — 400 cm -1 but with half the 
height corresponds to an antenna effect directly connecting sites \e\) 
and \e 3 ). Small deviations come from J23 being not zero. 



transferred out of the FMO complex [5, 6]. Experimental and 
theoretical studies have described two main relaxation path- 
ways in the FMO complex (HI [22l [23l EH). Here we shall 
focus on the dominant pathway that proceeds through site 1 to 
site 3, although our numerical analysis is based on simulations 
using the Modified Redfield approach as described in (22) of 
the entire eight- site FMO structure with an initial excitation 
on site 8. The sink is connected to site 3. Unless stated other- 
wise, the Hamiltonian parameters are taken from (27), and the 
spectral density of the environment is based on Adolphs and 
Renger fit to fluorescence narrowing experimental data (22) . 

Fig. [3^ shows the final population transferred to the sink 
for a simulation of a total time t = 2 ps at T = 77 K as a 
function of J\2 and £2. The transport efficiency shows a com- 
plex landscape of maxima and minima for the parameter re- 
gion explored, but we find that the Hamiltonian parameters 
of Ref. (27) (marked as a blue cross) are located at an op- 
timal point. Figure [3]) shows the energy level structure in 
both the local and exciton bases for the physical Hamilto- 
nian parameters of (27). We see in the exciton basis that the 
energy differences corresponding to the relaxation transitions 
connecting sites 1 , 2 and 3 are well matched to the maximum 
of the experimentally-determined spectral function, and there- 
fore this relaxation channel forms a striking example of an 
operational phonon antenna in a biological transport system. 
The inset of Fig. [3] shows p s [ n ^ for fixed Hamiltonian parame- 
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FIG. 3. a) Population of the sink, p s ir±(t — 2 ps) for different val- 
ues of the coupling between sites 1 and 2 (/12) an d the energy of 
site 2 (£2)- We clearly see that the actual values for J 12 and £2 lie 
at the top of a maximum. The inset shows /? s ink(/ — 2 ps) for the 
actual couplings as observed for the FMO, assuming an spectral den- 
sity function such as that of |29l but peaked at b) Local site 
and exciton energies for the H$ values corresponding to the FMO, as 
taken from |27|. The arrows show how the splitting might act as an 
antenna for the energy peak close to 180 cm -1 



ters and a varying peak position (Oh in the spectral density (as 
described in Appendix D). Ag ain, we find that p S mk is & multi- 
peaked function, but the most prominent maximum appears 
at ~ 190cm -1 , which is very close to the experimentally- 
fitted peak in the spectral function of Adolphs and Renger 
(180cm- 1 ) 1221 . 

The phonon antennae in a broader context — The phonon 
antennae principle which states that electronic interactions 
are tuned to achieve resonance with maxima of the environ- 
mental vibrational spectra and thus strongest response may 
be viewed in a broader context in which the tuning of elec- 
tronic degrees of freedom to match the frequency of external 
signals is exploited to achieve enhanced system response to 
facilitate sensing and transport. 

In fact, the phonon antennae principle is closely related to 
the Hartmann-Hahn condition l32ll which is of considerable 
importance in sensing devices. Here the challenge is that the 
sensor, e.g. an electron spin, may not be resonant with an ex- 
ternal system generating a signal (e.g. a different electron or 
nuclear spin). This problem may be overcome by continuous 
driving of the electron spin in the sensor at a strength that leads 
to a splitting of dressed states that matches the externally gen- 



erated signal and thus permits a strong response of the sensor 
(transitions between upper and lower dressed state) I33ti35ll . 

A further example in which tuning is not achieved by driv- 
ing is the recently proposed theory that olfaction (our sense 
of smell) is explained, at least partially, as originating from 
phonon assisted electron transport where electron source and 
drain in the receptor are separated by an energy splitting that 
is matched to a specific vibrational frequency of an odorant 
molecule. Electron tunneling, hence strong response, will oc- 
cur only if the molecular vibrations matches this energy gap 

£2Q. 

Discussion. — Taken together, these results suggest that the 
physically important relaxation pathway between sites 1 and 
3 is mediated by pigments which are spectrally and spatially 
positioned by the protein to efficiently sample the spectral 
function of the protein's fluctuations. Whether this optimal- 
ity is a determinant in the emergence of this structure in na- 
ture is beyond the scope of this study, but it is striking how 
well the phonon antenna concept can be used to rationalise 
the site energies and couplings of the pigments participating 
in this pathway. Other considerations, such as structural sta- 
bility and steric constraints imposed by the need to integrate 
the structure into the cellular environment are likely to drive 
the supramolecular structure of PPCs, however the multiplic- 
ity of local maxima we find does provide the opportunity to 
find an optimal working point for transport in the vicinity of a 
structure arising from these other factors. Indeed, we find that 
the FMO complex parameters correspond to a local maximum 
of the 1 — » 3 pathways efficiency. Interestingly, changes in 
temperature will also alter the inter-exciton energy that max- 
imises the relaxation rates in Eq. ([2]), and an evolving or adap- 
tive system might be able to change its coupling structure to 
maintain the inter-exciton energy differences at the maximum 
of the thermal noise strength. 

We also point out that while the Redfield formalism al- 
lows us to give a physically clear description of the essential 
phonon antenna mechanism, the use of Markovian Redfield 
theory for PPC problems neglects a number of important dy- 
namical effects (environmental memory, stokes shift, discrete 
modes in the spectral function) which are crucial for under- 
standing the persistence of inter-eigenstate coherence that is 
generated in experiments fT6ifT8ll . These effects have been 
neglected here for simplicity, but exploring these effects in 
the context of the phonon antenna mechanism could lead to an 
even more sophisticated principle, involving transient, broad- 
band sampling of the spectral function and evolving exciton 
spectra and localisation lengths which could, potentially, ac- 
tivate or deactivate relaxation pathways dynamically. Finally, 
the concept of using coherent interactions to alter the sampling 
of a (fixed) environmental spectral function could be readily 
achieved in a variety of atomic and condensed matter systems 
El El, and could also be used as a sensitive probe for mea- 
suring spectral functions by measuring transport rates as func- 
tions of a controllable coupling parameter. 
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SUPPLEMENTARY INFORMATION 

Appendix A: Detailed Model 

The total Hamiltonian (fi = 1) describing electron-phonon 
interactions consists of the system Hamiltonian Hg which in- 
cludes the electronic degrees of freedom, the free Hamiltonian 
Hr of the phonon bath describing fluctuations in the protein 
structure, considered as a set of N independent bosonic baths 
{ / = 1 , 2. . .N} coupled to each of the sites (chromophores), and 
the linear coupling of fluctuations to the excitations V : 

H=H S + H R +V, H S = J> + J^Jij 

Hr =Y*ll COka ik a iki V = + lO^I ( 4 ) 

i k ik 

where a^a\ k are the usual creation and annihilation operators, 
respectively, satisfying [a^^a]^ = 8^. For simplicity, we will 
assume an identical coupling to the environment for each site 
gik = gk with /{(to) = L*gjfc5(fi) - co k ) being the vibrational 
spectral density. 

In the exciton basis, Hg is diagonal, Hs\e n ) = E n \e n ), 
and we will write |/) = £ n Q \e n ). In that basis the Hamil- 
tonian can be expressed as H$ = Y. n ^n\e n ) {e n \ and V = 



Hik8ikC l m C l *(a ik ^aJ k ) \e m ) (e n \, where we have ne- 
glected the Lamb and Stark shifts corresponding to m = n in 
the last term. 

A state |0), or sink, has been introduce to model the trans- 
fer to the RC, where the excitation will end up eventually 
as a relaxation process coming from site 3, so the behaviour 
of the system will account for an additional channel l^^o = 
1/(1 ps) . Although not fully equivalent, when using the mas- 
ter equation in the exciton picture we will model this irre- 
versible transfer to the sink as happening at rate T^^o = 

|(^|3)| 2 r 3 ^ . 

In all cases we study the time evolution of the sink's popu- 
lation following the selective excitation of site a 9 \\j/) = \a). 
In the considered 3 sites model system, a = 1, while in the 8 
site model of FMO, it is a = 8. 



Appendix B: Non-Markovian analysis of a three site network 

In this section we analyse the 3 -site case beyond the con- 
cerns exposed in the markovian treatment, with the aim to 
show how manifestations of the phonon antenna also occur in 
a paradigmatic example of non-markovianity. For that matter 
we will assume that all levels in the network are strongly cou- 
pled to identical independent harmonic oscillators, all tuned to 
the same frequency (Oh, representing therefore the coupling to 
a unique discrete mode. The coupling between sites would be 
given by the same matrix Hg. 

This model presents the advantage of enabling a direct anal- 
ysis in the local site basis, which lets us avoid making the ap- 




FIG. 4. Evaluation of non-Markovian effects. Local sites are coupled 
to damped harmonic oscillators, so that the reduced dynamics for the 
electronic degrees of freedom alone is in general non reproducible 
with a Markovian master equation. 



proximation T en ^$ = | (e n |3) | 2 T3^o for the transfer rate to the 
sink node. The local Hamiltonian described in Eq. ^ fits this 
model if we assume that the sum of &'s restricts to just one 
mode. The dissipation of each mode will be modelled using 
a master equation formalism, where the harmonic oscillators 
relax at the same rate (r = 60 cm -1 ) to a (Markovian) thermal 
bath at temperature T, according to the Lindblad superopera- 
tor: 



=2diss = J^r[n((0 H ) + 1] [ciipa] - ^ajdip - ^pajcii] 

+ J^T[n((OH)][aJpai - ^a t ajp - ^pata]} (5) 

The coupling constant g is chosen to match the quasi 
Lorentzian profile previously presented in l26l . There, it is 
shown that a two level system coupled to a bath with that pro- 
file behaves equivalently as if it were coupled to a damped 

oscillator with coupling constant g = (Oh \[^ — \ V^(Oh ~ 
80 cm -1 Within this approach, the oscillator is damped within 



an environment of Ohmic nature ^f((o) = ^0). Although 
that is not exactly our case, it can be used as a first approxi- 
mation to study the impact of non-Markovian features. 

The transfer to the sink is modelled as an irreversible pro- 
cess from site 3 to site with T3^o = 1 / ( 1 P s ) as described 
by: 



-^sink — T3->0 



|0)(3|p|3)(0|-i|3)(3|p-ip|3)(3| 



(6) 



The inset of Fig. [5]shows some changes in the phenomenon. 
It takes place close to the expected value, although the shift is 
prominent. The explanation relies on the direct coupling to the 
harmonic oscillators, which is so strong that leads to a mod- 
ification in the excitonic levels of the whole system, the final 
configuration being way more complicated than the one ex- 
plained in section II. The concept remains valid, but of course 
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FIG. 5. Simulation of 10 ps transfer with oscillators tuned at 200 
cm - 1 at T=4 K. The actual couplings lie close to a maximum. The 
inset shows how a phonon antenna effect appears if the mode is tuned 
at roughly 237 cm" 
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Appendix D: Quantitative Figure of merit 

We have envisioned a figure of merit to check how well can 
the concept of phonon antenna alone predict the main features 
of the landscape as shown in Fig. [3^. For the same variations 
of J\2 and £2 we have studied the level structure and computed 
the simple value: 



the intuition about the configuration of the antenna is not so 
precise. 

Fig. [5] shows a plot with the mode fixed at 200 cm -1 . As 
we see, the values of the Hamiltonian are close to a maxi- 
mum in transfer yield showing once more the forming of an 
antenna. 



Appendix C: Spectral function 

For the sake of modelling appropriately the spectral func- 
tion used for the FMO case, we have relied on the function 
proposed by Adolphs and Renger l22l . In that reference a 
combination of super-ohmic densities is used and then a sin- 
gle effective high-energy vibrational mode H is added with 
Huang Rhys factor S H = 0.22 and co H = 180 cm -1 . In our 
case, instead of using a delta function, we introduce a pure 
Lorentzian with the same characteristics but a linewidth of 
Fh = 60 cm -1 . With those prescriptions the spectral function 
can be written as 



^antenna = max { - | (0 H - \E+ - £_ 1 1 - \C0 H - \E--E G \\, 

- l -\2(o H -\E+-E G \\} (8) 

The higher the value of Antenna, the more assisted by the envi- 
ronment the transition to the sink would be, as inferred from 
the scheme described. Fig. [6] shows the map of Antenna m the 
same range of parameter values as Fig. [3] 
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FIG. 6. F an tenna for varying values of Ju and £2 with coh = 180 cm 
as in Fig. [3^ 



